The glass transition in a model silica glass: 
evolution of the local structure 



P. Jund and R. JuUien 

Laboratoire des Verres - Universite Montpellier 2 
Place E. Bataillon Case 069, 34095 Montpellier France 

Abstract 

We use molecular dynamics simulations and the Voronoi' tessellation to study 
the geometrical modifications as a function of temperature in a model silica 
glass. The standard deviation of the cell volumes, which is a measure of 
the local density fluctuations, decreases with decreasing temperature, as if it 
would like to vanish at zero temperature. This evolution towards an ordered 
state is frozen out at the glass transition and consequently an amorphous 
sample is obtained at low temperature. This structural freezing following 
upon the glass transition is noticeable in all the other geometric characteristics 
of the Voronoi cells and a possible interpretation in terms of geometrical 
frustration is proposed. 
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I. INTRODUCTION 



Silica is a common material which is of great importance in chemistry, geology and 
industrial applications. It is also a prototype of a network forming glass. All these reasons 
explain why it has been the topic of a great amount of studies. Nevertheless many features 
of this typical "strong" glass need still to find a satisfactory explanation. For example the 
origin of the First Sharp Diffraction Peak (FSDP) is still controversial (EUiott 1991; Gaskell 
and Wallis 1996; Fayos, Bermejo, Dawidowski, Fischer and Gonzalez 1996). The origin 
(and the connection to the FSDP) of the so-called Boson peak remains the topic of many 
studies, both experimental and theoretical (Shuker and Gammon 1970, Martin and Brenig 
1974; Akkermans and Maynard 1985; Buchenau, Zhou, Niicker, Gilroy and Phillips 1988; 
Mahnovsky, Novikov, Parshin, Sokolov and Zemlyanov 1990; Novikov and Sokolov 1991; 
Sokolov, Kishuk, Soltwisch and Quitmann 1992; Borjcsson, Hassan, Swenson, Torell and 
Fontana 1993; Gurevich, Parshin, Pelous and Schober 1993; Achibat, Boukenter and Duval 
1999; Schirmacher and Wagener 1993; Bermejo, Criado and Martinez 1994; Terki, Levelut, 
Boissier and Pelous 1996). With the development of the computing speed a new type of 
studies has emerged in the past decade which can be called "computer experiments" . Indeed 
numerical simulations have been used to study the vibrational spectrum (Jin, Vashista, 
Kalia and Rino 1993; Taraskin and Elliot 1997a; Horbach, Kob and Binder 1997; Taraskin 
and Elliott 1997b; Guillot and Guissani 1997) or the structural characteristics (Woodcock, 
Angell and Cheeseman 1976; Holender and Morgan 1991; Delia Valle and Andersen 1992; 
Servalli and Colombo 1993; Nakano, Bi, Kalia and Vashista 1993; Sarnthein, Pasquarello 
and Car 1995; VoUmayr, Kob and Binder 1996) of various model silica glasses. Within this 
framework we present in this paper a classical molecular dynamics study essentially focused 
on the evolution of the local structure in such a system. 

The first step in that kind of endeavor is to choose the interparticle potential: we decided 
to use the two-body potential proposed by van Beest et al. (van Beest, Kramer and van 
Santen 1990). Indeed a study of the influence of the quenching rate on the properties of 
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amorphous silica has shown that this potential based on ab initio parameters gives excellent 
results compared to experimental data both for the structural characteristics (VoUmayr et al. 
1996) and the vibrational properties (Taraskin and Elliott 1997b) of vitreous silica, even if, 
in simulations the quenching rates are always faster than in real systems . Similarly to what 
has been done earlier in the case of high pressure silica samples (Rustad, Yuen and Spera 
1991) we combine the molecular dynamics with the Voronoi tessellation scheme in order 
to have a better insight into the structural evolution of the system during the quenching 
procedure. We want to know if and how the glass transition leaves a "signature" in the 
geometrical characteristics of the Voronoi cells. We recall that the Voronoi cell attached 
to a particle is an extension, for disordered systems, of the Wigner-Seitz cell and gives 
information on the local structure around this particle. 

Our results as a function of temperature show the structural freezing consecutive to the 
glass transition in all the geometric characteristics of the Voronoi cells. We suggest that the 
structural freezing below Tg is mainly a consequence of the impossibility to build an ideal 
structure as a consequence of the so-called "geometrical frustration" (Kleman and Sadoc 
1979; Kleman 1983, 1985, 1989; Sadoc and Mosseri 1997). 

In section II a description of the simulations will be given. We present our results in section 
III and discuss them in section IV. Finally in section V we present our conclusions. 

II. SIMULATIONS 

If one wants to perform realistic simulations on silica glasses the major point is to choose 
an interaction potential which gives reasonable results compared to experimental data. Sev- 
eral choices are today possible, but in the last years one of the most successful classical 
potentials is the so-called BKS potential developed by van Beest et al. (van Beest et al. 
1990). Though designed originally from the crystalline phases of silica, it has been shown 
that it also describes very well the structural (VoUmayr et al. 1996) and vibrational (Taraskin 
and Elliott 1997b) properties of amorphous sihca. 
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The functional form of the pairwise BKS interaction between two particles i and j is given 



where rij is the interparticle distance, e the charge of an electron and the parameters Aij, 
Bij and C^j are fixed as follows: Agio = 18003.7572 and Aqo = 1388.773 eV]Bsio = 4.87318 
and Boo = 2.76 A"^; Csio = 133.5381 and Coo = 175.0 eVA^. Note that except for the 
Coulomb interaction {qsi — 2.4 and qo — —1-2) there is no interaction between Si atoms. 
This original form contains an unphysical property at short distances since it diverges to 
minus infinity. To overcome this drawback which is especially annoying at high temperature, 
we have added a short range repulsive term (~ l/r^°) which insures that the potential goes 
to infinity at small interatomic distances and is practically equivalent to the original poten- 
tial for Tij > 1.2A (Si-0 interaction), and Vij > 1.6A (0-0 interaction) without introducing 
any artificial energy barrier. 

The Coulomb interactions were computed using the Ewald summation method (Ewald 1921; 
Allen and Tildesley 1990) with a characteristic constant k, — 6.5/L, where L is the cubic 
box size, and considering 501 k- vectors in reciprocal space < 6 x 27r/L). These values 
insure that the potential energy is obtained with a relative error smaller than 5 x 10~"^%. 
No cut-off was used for the pairwise interaction. 

We performed molecular dynamics simulations for microcanonical systems containing 216 
silicon and 432 oxygen atoms confined in a cubic box of edge length L = 21.48A, which cor- 
responds to a mass density of ~ 2.18g/cm^ very close to the experimental value of 2.2g/cm^. 
Periodic boundary conditions were used to limit surface effects. In order to insure energy 
conservation even at high temperature a timestep of 0.7 fs was necessary. This value is sub- 
stantially lower than the one used in previous studies (1.6 fs (VoUmayr et al. 1996) or 1.0 fs 
(Taraskin and Elliott 1997b)) which may be due to our "conservative" potential correction. 
The 4-th order Runge-Kutta algorithm was used to integrate the equations of motion. 
The glass configurations were obtained by quenching well equilibrated initial liquid samples 
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by 




(1) 



obtained by melting /9-cristobalite crystals at a temperature around 7000 K. After full equi- 
libration of the liquid 40000 timesteps), the system was cooled to zero temperature at a 
quench rate of 2.3 x 10^"^ K/s which was obtained by removing the corresponding amount 
of energy from the total energy of the system at each iteration. Due to computer time 
limitations this cooling rate is rather fast but it has the advantage compared to other pro- 
cedures involving either stepwise cooling (VoUmayr et al. 1996) or temperature dependent 
rates (Taraskin and EUiott 1997b), to be hnear and continuous all along the quenching pro- 
cedure. At several temperatures during the quenching process the configurations (positions 
and velocities) were saved. Each configuration was used to start a constant-energy molecular 
dynamics calculation during which the temperature was recorded as a function of time. The 
temperature was in all cases remarkably stable and only slight relaxation effects could be 
observed. Nevertheless to avoid transient configurations we allowed for each temperature 
and for each sample, 10000 relaxation steps followed by 50000 supplemental time steps (for 
a total simulation time of 42 ps) during which all the calculations were done. Apart from 
the calculation of standard quantities (radial pair distribution function, mean square dis- 
placement) we included also in our molecular dynamics code a Voronoi' tessellation scheme 
similar to the one that we have developed for monocomponent soft-sphere glasses (Jund, 
Caprion and JuUien 1997a; 1997b). This scheme has been modified to take into account sev- 
eral types of atoms and thus it permits to follow the local structure around the silicon atoms 
and the oxygen atoms as a function of temperature during the quenching procedure. Here 
the Voronoi cell is always defined as being the region of space closer to a given atom center 
than to any other and no dissymetry between the two components has been introduced as it 
is the case in the "navigation map" procedure (Gellatly and Finney 1982). The Voronoi cell 
characteristics as well as all the other quantities have been averaged over samples obtained 
from 5 independent starting configurations in order to improve the statistics of the results. 
The whole simulation lasted for more than 3 million timesteps which were run on 4 nodes 
of an IBM/SP2 parallel computer. 
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III. RESULTS 



As said earlier this potential has already been used in other studies of amorphous silica 
and since we did not modify the BKS parameters our results concerning the radial pair 
distribution (VoUmayr et al. 1996; Taraskin and Elliott 1997b) or the diffusion constant 
(Guillot and Guissani 1997) are exactly identical to the referenced results and therefore we 
do not come back to these standard results here. Our aim is to locahze the glass transition 
temperature Tg through the study of the structural characteristics (via the Voronoi tessella- 
tion) of our model silica system. Nevertheless a straightforward way of determining Tg is to 
monitor the potential energy versus the temperature, as has already been done for another 
model silica glass (Delia Valle and Andersen 1992). The evolution of the average potential 
energy per particle versus the temperature is shown in Fig. 1. With increasing temperature, 
the potential energy as well as the standard deviation increases as expected. Nevertheless 
one can observe an acceleration of this increase between 3000 and 4000 K, corresponding to 
the passage from a sohd to a liquid behavior. Due to the fast cooling rate the value of Tg 
is much higher than the experimental value (1446 K (Angell 1988)) but it is coherent with 
the value of ~ 3500 K obtained from the fit of Tg versus the quenching rate proposed by 
VoUmayr et al. (VoUmayr et al. 1996). Note that we observe a non negligible increase of 
the potential energy in the glass phase contrarily to what was obtained by Delia Valle et 
al (Delia Valle and Andersen 1992, Fig. 3) who used a direct minimization procedure after 
the quench to investigate the "inherent structures" (StiUinger and Weber 1982). Here we 
let the system evolve freely after the quench and it seems that even at low temperatures 
structural relaxation occurs. Note also that for such a fast cooling rate the structure of the 
resulting silica might very well be slightly different than the one observed experimentally, 
since it is known that to get the most "ideal" glass, one should use the slowest quenching rate 
which does not produce crystallization (Kauzman 1948; Jund, Caprion and JuUien 1997c). 
However we recall that previous authors (VoUmayr et al. 1996; Taraskin and Elliott 1997b), 
using the same potential with a quenching rate of the same order of magnitude as ours, were 
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able to reproduce several physical properties of real silica. 

Once we know approximately the value of Tg (to our purpose this level of accuracy is suf- 
ficient) we can tackle the study of the geometrical characteristics of the Voronoi cells in 
order to follow the local structure as a function of temperature. All the characteristics of 
the cells have been obtained (surface, number of faces, number of edges, etc..) but we want 
to discuss here only some representative quantities. The first one is the variation of the 
volume of the Voronoi cells. This variation is represented in Fig.2 for the silicon (a) and the 
oxygen (b) atoms together with the corresponding standard deviations (c). With decreasing 
temperature the volume of the silicon cell decreases while the volume of the oxygen cell 
increases (these opposite variations are a consequence of our constant- volume calculations) . 
Again a change of behavior is visible and corresponds to a slowing down of the evolution 
below the glass transition temperature. An even more striking behavior is observed in Fig. 2c 
where the standard deviation ay is plotted as a function of temperature. This is a quantity 
of physical interest since it measures the local density fluctuations around the particles. For 
both types of atoms, ay decreases with decreasing temperature as if it would like to tend 
to zero at T = OK but then below 4000 K this trend is stopped and finally uy saturates 
around non-zero values characteristic of spatial disorder. This is a direct observation of 
the low-temperature saturation of the density fluctuations which is a signature of the glass 
transition. Even though not reported here, exactly the same behavior can be observed for 
the evolution of the surface of the Voronoi cells as a function of temperature. 
To investigate further the structural evolution during the quench, wc looked at the angle dis- 
tributions. To numerically determine the bond angle 0-Si-O (resp. Si-O-Si), we determine 
for each Si (resp. O) atom the two nearest O (resp. Si) atoms and we calculate the angle 
between the two corresponding segments Si-0 (resp. 0-Si), the result being averaged over 
all the Si (rep. O) atoms. First we studied the tetrahedral 0-Si-O angle which should be 
ideally equal to 109°. 47 in a perfect tetrahedron. As can be seen in Fig. 3a, this angle varies 
between 110°. 5 at low temperature and 117° at 7000 K, with a slight change of behavior 
around Tg. The same behavior is observed for the corresponding standard deviation (Fig.3c 
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full circles) which is quite small and decreases when T decreases. This shows firstly that 
with increasing temperature the Si04 tetrahedra survive even in the liquid phase but become 
more and more distorted and secondly that the glass transition does not strongly affect the 
local environment around the silicon atoms. On the contrary the glass transition is more 
clearly visible in the angle Si-O-Si, which measures the relative position and orientation of 
two neighboring Si04 tetrahedra, as can be seen in Fig. 3b. With decreasing temperature 
within the liquid phase, the Si-O-Si angle increases and seems to converge towards 180° (a 
least square quadratic fit of the four points in the liquid phase gives an extrapolated value 
of 175°), but again below 4000 K this increase slows down and finally the angle converges 
towards a value close to 150°, a value coherent with previous simulations (VoUmayr et al. 
1996; Taraskin and Elliott 1997b), but slightly higher than the value 144° found in X-ray 
diffraction experiments (Neuefeind and Liss 1996). The decrease of this angle is coherent 
with the views of a densifying network with increasing temperature this densification taking 
place around the oxygen atoms. Also since it measures the relative orientation between 
two neighboring tetrahedra, this decrease corresponds to a decrease of the effective volume 
of the oxygen atoms with increasing temperature and since we work at constant volume it 
implies an expansion of the silicon volume, which is indeed the behavior observed in Fig. 2. 
Concerning the standard deviation of the Si-O-Si angle represented in Fig. 3c (open circles) 
it increases rapidly in the liquid phase (above 4000 K) while in the glass phase this increase 
is more slow. This again is an illustration of the structural freezing below Tg. 
Another quantity which reflects the glass transition is the effective coordination number, z, 
which is in fact the average number of faces of the Voronoi cells. This quantity should be 
considered with some care since, generally, it does not correspond to the true "chemical" 
coordination number. To give an example, in /3-cristobalite the Voronoi cell of the silicon 
atoms is a tetrahedron [z — 4) while the Voronoi' cell of the oxygen atoms is the polyhedron 
represented in Fig. 4 with z — d,. In the Voronoi "sense" the nearest neighbors of an oxygen 
atom in /3-cristobalite arc two silicon atoms (represented by the triangular faces in Fig. 4), 
and six oxygen neighbors (represented by the pentagonal faces in Fig. 4). Such a cell is 
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a truncated double-tetrahedron and is, in fact, highly "degenerate": some of its edges are 
shared by four neighboring cells (instead of three usually) . An infinitely small random per- 
turbation of the atomic positions (~ 10~^A) is sufficient to create new extra-small triangular 
faces and the coordination number of the oxygen atoms flips to 2; ~ 19.7 while the coordi- 
nation of the silicon atoms remains close to 4. This could be avoided by doing a smoothed 
Voronoi tessellation which excludes small faces with an area smaller than 10% of the largest 
cell face (Rustad et al. 1991) but we have rather chosen not to use this technique and ana- 
lyze raw data. The variation of 2; as a function of temperature is represented in Fig. 5a for 
the silicon atoms and Fig. 5b for the oxygen atoms. The glass transition can again clearly 
be identified: for the silicon atoms, z increases with increasing temperature with a slope 
that is more important above 4000 K. For the oxygen atoms, z is approximately constant in 
the glass phase and then decreases for temperatures above Tg. It should be noted that the 
relative variation of z between and 7000 K is small, especially for the oxygen atoms. In 
the liquid phase, with decreasing temperature, the local structure around the silicon atoms 
evolves towards a perfect tetrahedral arrangement [z = 4) and then converges to a value 
higher than 4 in the glass phase due again to the structural freezing below Tg. 



IV. DISCUSSION 

All these results can be discussed in the light of the conclusions drawn for the analo- 
gous geometrical analysis performed in a monoatomic soft-sphere glass (Jund et al. 1997a, 
1997b). In that case, the system tries to reach an icosahedral arrangement (with dodecahe- 
dral Voronoi cells) when the temperature is lowered down from the liquid phase, but, since 
such an arrangement can not be realized at large distances in the regular three-dimensional 
space, the system gets frozen in a glass phase below a characteristic glass temperature. Such 
geometrical frustration effects were in that case the consequence of the degeneracy between 
the face-centered-cubic (FCC) and the hexagonal-close-packed (HCP) structures (Sadoc and 
Mosseri 1997). 
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Since all the standard deviations presented here seem to go to zero with decreasing tem- 
perature, when extrapolated from the liquid phase, it is reasonable to assume that in the 
case of silica also a T = unreachable ideal local structure exists. It is also reasonable to 
assume that such an ideal arrangement corresponds to a perfect tetrahedral order for the 
four oxygens bounded to a given silicon atom, as it is for almost all of the known crystalline 
structures of silica. This assumption is supported by the behavior of z, which tries to ex- 
trapolate to 4 when T is lowered from the liquid phase, and it is not incompatible with our 
results for the variation of the angle 0-Si-O with temperature (see Fig. 3a). Even if this angle 
does not show a major change of behavior at T^, an extrapolated value of 109. °47 at T = 
is not inconsistent with the reported data above Tg. Moreover, since the angle Si-O-Si seems 
to extrapolate to 180°, one can imagine that the ideal structure is made of tetrahedral units, 
like the sp^ coordination of carbon where the silicon atoms would be located at the carbon 
places and oxygen atoms located in the middle of the C-C bonds (see Fig. 6a). Indeed the 
tendency to build such a local structure should result from the form of the potential. In 
particular the tetrahedral arrangement of the oxygens around a silicon atom results from a 
combination of the Si-0 attraction and the repulsion between the oxygens. The tendency 
to align the Si-O-Si bridges between neighboring tetrahedra is more subtle however, since 
the long range nature of the ionic part of the potential certainly plays a role. 
It is interesting to notice that, among all the known crystalline structures of silica, two 
particular structures (at least) fully satisfy these criteria, namely the /3-cristobalite and the 
tridymite structures. In these two structures the above defined tetrahedral units arc stacked 
with sequences ABCABC...(Fig.6b) and ABABAB...(Fig.6c), respectively, hke in FCC and 
HCP structures. In fact these structures can be simply built from FCC and HCP struc- 
tures, by adding to the original structure another one shifted by a fourth of the diagonal 
of the cubic cell (in the FCC case) and by 3/8 of the c-axis of the hexagonal cell (in the 
HCP case). They correspond respectively to the diamond and wurtzite structures of carbon. 
When these two structures arc considered with the same density, they have exactly the same 
Si-0 distance dsiOi and therefore, in the two cases the Voronoi cells for the silicon atoms are 
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regular tetrahedra with the same volume Vsi — V^dg^Q. The volume of the oxygen cells, Vq, 
is also the same and therefore the two structures are characterized by the same Si/0 volume 
ratio R = Vsi/{2Vo) = 9/55 = 0.164, independent on the density. In Fig.7 we have plotted 
R — Vsi/{2Vo) as a function of T, as calculated from our simulations, and reported the 
value R — 0.164 at T = (open circle). When decreasing the temperature from the liquid 
phase R decreases as if it would like to reach a value quite close to (or even lower than) 
0.164. Therefore, one could conclude that may be the ideal structure that the system would 
try to realize when lowering temperature would be a mixture of tridymite and cristobalite 
phases. In fact this is not a satisfying solution since these two phases can very well arrange 
themselves under the form of crystallites separated by dislocations or stacking faults, and 
therefore crystallization should actually occur in our system. 

One can find a solution by pursuing the analogy with soft sphere glasses or hard sphere 
packings for which also two degenerate ground state arrangements FCC and HCP exist. In 
that case the system gets "frustrated" by trying to reach a local icosahedral structure (Jund 
et al. 1997a, 1997b; JuUien, Jund, Caprion and Quitmann 1996) which does not exist in the 
flat regular space. Similarly to these model systems, the frustration can be resolved by con- 
sidering in our case a curved space with positive curvature, namely the sphere 5*3 (Sadoc and 
Mosseri 1997). Consider a regular tetrahedron with a sihcon atom in the middle and oxygen 
atoms at the centers of the faces. This tetrahedral unit contains one sihca molecule. One can 
exactly tile an 5*3 space with 600 tetrahedra like this, and the resulting Si02 structure (which 
contains 600 silicons and 1200 oxygens) satisfies all the local requirements defined above. 
The 120 vertices of the unit tetrahedra are located on the so-called {3,3,5} polytope (Sadoc 
and Mosseri 1997). Note that an important difference with tridymite and cristobalite is that 
this ideal structure is made of five-membered ..Si-0-... rings, instead of six-membered ones, 
since there are five unit tetrahedra around a tetrahedron edge in the {3,3,5} polytope. It is 
not trivial to calculate the exact value of the ratio R for this structure but a good approx- 
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imation can be obtained by considering the tetrahedron unit in the flat three dimensional 
space tangent to the curved space at the tetrahedron center. The sihcon Voronoi' cell, limited 
by the bisector planes of the Si-0 bonds, is a tetrahedron of volume 1/8 of the volume of 
the unit. Therefore the remaining volume for the two oxygen atoms is 7/8 and consequently 
R— 1/7 = 0.143. This value is only approximate due to the fact that volumes do not exactly 
scale as distance-cubes in the curved space 5*3 of positive radius r^. An estimation of the 
error made can be obtained by performing the same reasoning for two concentric spheres of 
radii r and r/2, for which one knows analytic formulae for their volumes (Sadoc and Mosseri 
1997). Performing an expansion in r/r^, one gets R ~ (1/7)(1 + (2/35)(r/rs)^...) instead of 
1/7. Taking for r the radius of the sphere of same volume than the one of the unit tetrahe- 
dron (which is the volume of 5*3 divided by 600), this formula gives 0.144 instead of 0.143. 
Hence, in view of the roughness of the extrapolations done, such an error can reasonably be 
neglected. 

This value is represented by the open square in Fig. 7. It appears that this value corre- 
sponds to a better extrapolation of the four points above 4000K than the one obtained from 
tridymite or cristobalite (the second order flt of the four liquid points leads to a value of 
0.136). Therefore one is tempted to conclude as in the soft sphere case (Jund et al. 1997a; 
1997b): when lowering the temperature from the liquid phase the system evolves as if it 
would try to build locally such an ideal structure, but since this structure cannot be real- 
ized in the regular three dimensional space, the systems gets frozen in a glassy state below 
Tg. Obviously, the above interpretation should be considered as a suggestion and there are 
certainly other ways to resolve the frustration. One way is to consider the "polytope 240" 
described in chapter 2 of the Sadoc and Mosseri book (Sadoc and Mosseri 1997). In that 
case, one would get six-membered rings. One could also imagine solutions with seven mem- 
bered rings. The only argument in favor of our simple suggestion is the nice extrapolation of 
the volume ratio data to 0.143. We have checked that this ratio is signiflcantly larger for the 
polytope 240 solution. Definitely more calculations are needed to validate our proposal. In 
particular a systematic enumeration of the n-membered rings as a function of temperature 
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would be helpful. Such an enumeration already exists in the literature (VoUmayr et al. 1996) 
but has been done as a function of quenching rate. 

V. CONCLUSIONS 

With the use of classical molecular dynamics simulations combined with the Voronoi tes- 
sellation we have studied the evolution of the local structure around the particles of a model 
silica glass as a function of temperature. The glass transition temperature, T^, obtained 
from the variation of the potential energy versus the temperature is coherent with the value 
expected for systems quenched at 2.3 x 10^^ K/s. This transition is clearly visible in the 
evolution of all the geometric characteristics of the Voronoi cells. The study of the local den- 
sity fluctuations clearly demonstrates the saturation of these fluctuations below Tg which 
is the usual signature of the glass transition. The study of the tetrahedral 0-Si-O angle 
shows that with increasing temperature the Si04 tetrahedra survive but become more and 
more distorted. The variation of the angle Si-O-Si between two corner-sharing tetrahedra 
is coherent with the views of a densifying network with increasing temperature, this densi- 
fication happening around the oxygen atoms, while a volume expansion occurs around the 
silicon particles since we are considering microcanonical ensembles. The different geometric 
characteristics and in particular the average coordination number evolve with decreasing 
temperature as if the system would like to reach an ideal structure at T = 0, which cannot 
be realized due to geometrical frustration. This evolution is frozen below the glass transition 
temperature, and finally the system converges towards an amorphous structure. 
Our results show that the use of the Voronoi cell characteristics gives not only useful in- 
formations on the local structure but can also be used to determine the glass transition 
unambiguously. This means that even if nothing dramatic happens at the glass transition 
concerning the local structure, something happens, which is basically a dynamic freezing of 
the natural evolution of the structure towards an unreachable ideal structure. Of course we 
are far away from the timescales used in experiment, but this study, together with others. 
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permits to investigate what happens on the microscopic level in an attempt to explain what 
is observed at the macroscopic level. 
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FIGURES 



Fig. 1. Average potential energy per particle as a function of temperature. 

Fig. 2. Variation of the cell volume versus temperature: (a) silicon; (b) oxygen, 
(c) Standard deviation ay versus temperature: •: silicon; o: oxygen. 

Fig. 3. Variation as a function of temperature of: (a) 9, the 0-Si-O angle, (b) 0, the Si- 
0-Si angle; (c) the standard deviations and cr^ versus temperature: •: silicon; o: oxygen. 

Fig. 4. The Voronoi cell of an oxygen atom in /3-cristobalite. 

Fig. 5. Variation of the effective coordination, z, as a function of temperature: (a) silicon 
atoms; (b) oxygen atoms. 

Fig. 6. The ideal tetrahedral unit around a silicon atom (a) is shown together with 
the trydimite (b) and /3-cristobahte (c) structures viewed from the top. In (b) and (c) two 

successive layers connected to a single top unit are represented (the positions of the oxygen 
atoms have been omitted). 

Fig. 7. Variation of R = Vsi/2Vo as a function of temperature, o: value of R obtained 
in tridymite or /9-cristobalite; □: value of R obtained in the ideal structure built on S3. 
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